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. As is well known, energy cost can greatly impact the deployment of battery-powered sensor networks 

in remote environments such as rivers or oceans. Motivated by this, we propose here an energy-based 
metric and associate energy-based Voronoi partitions with mobile vehicles in constant flows. The metric 
corresponds to the minimum energy that a vehicle requires to move from one point to another in the 
I flow environment, and the resulting partition can be used by the vehicles in cooperative control tasks 

(~| , such as task assignment and coverage. Based on disk-based and asymptote-based approximations of the 

, Voronoi regions, we determine a subset (or lower bound) and superset (or upper bound) of an agent's 

Voronoi neighbors. We then show that, via simulations, the upper bound is tight and its cardinality 
remains bounded as the number of generators increases. Finally, we propose efficient algorithms to 
compute the upper bound (especially when the generators dynamically change), which enables the fast 
' calculation of Voronoi regions. 

o ' 
p , 

I. Introduction 

Due to the proliferation of low-cost sensing, communication, and computation devices, large 
groups of mobile vehicles equipped with sensors can be deployed into flow environments (e.g., 
>< : rivers, lakes, oceans) to efficiently perform monitoring tasks. Depending on the nature of the 
application, multiple mobile vehicles can be coordinated based on different objectives. For search 
and rescue missions, a priority is to find/reach the target within the shortest time. However, for 
non-urgent tasks such as the monitoring of harmful algae blooms, maximizing the lifetime of the 
whole group of mobile vehicles can be more critical, as mobile vehicles are commonly powered 
by batteries with limited capacity. This motivates the study of minimum energy cooperative 
control algorithms for mobile vehicles in flow environments. 

In this paper, we study a Voronoi partition associated with the minimum energy required for 
a vehicle to move from one point to another in a constant flow environment. We first derive 
an explicit expression for the energy-based metric, and study the Voronoi partition based on 
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this metric using the vehicle locations as the set of generators. Similar to the time metric 
counterpart Jil, the Voronoi partition can then be used in the design of efficient target-assignment 
(or task allocation) algorithms (e.g., some other work |[2l, Q). By assigning a vehicle to the 
targets that fall into its Voronoi region and guiding its motion appropriately, the vehicles can 
minimize the average energy spent by the group in servicing stochastic tasks that arrive according 
to a slow-rate Poisson distribution. However, contrary to the Euclidean case [IJ, the Voronoi 
region defined by a general metric can be very involved. On the other hand, upper and lower 
approximations of the regions can be just enough to implement a coverage or target assignment 
algorithm; see IIH. Motivated by this, we propose methods to bound the set of Voronoi neighbors 
of a vehicle, which simplifies the calculation of Voronoi cells by vehicles. These are based on 
the following considerations: (i) the characterization of Voronoi region boundaries as hyperbolas, 
(ii) approximations of Voronoi regions by means of circles and polygons, and (iii) the derivation 
of a simple test (refer to Theorem [2l Corollary [U and Proposition ^ that allows to discard 
vehicles that cannot be Voronoi neighbors. The test leads to an upper bound on the set of 
Voronoi neighbors of a vehicle. By generating vehicle locations (i.e., the set of generators 
for the Voronoi partition) independently according to a uniform distribution, we show that the 
average number of generators in the upper bound is bounded (via simulations) by 4.5. Since 
the set of generators in the upper bound is sufficient for calculating Voronoi cells, the approach 
based on this upper bound (instead of using all generators) can save significant amount of time 
when calculating Voronoi cells, especially for applications with large amount of mobile vehicles. 
Therefore, we propose different algorithms to calculate the upper bound, especially when dealing 
with constantly changing generators. 

Previously, Voronoi partitions in flow environments have been studied in connection to the 
shortest traveling time metric Q-lIll. In contrast, there are relatively fewer works on the energy 
metric [[T0ll - [[T4ll . For example, in ifTTl - lfTBl . the goal is to find a path with the minimum energy 
loss between given source and destination points in piecewise constant regions; in [[T4| , the 
minimum energy metric (which is the same as this work) is used but the flow is modeled 
as a quadratic function (in this case, there is no explicit expression for the energy metric). 
In terms of approximating Voronoi cells, the work in dU, [[T5l propose methods to deal with 
Voronoi partitions induced by the Euclidean distance metric (called standard Voronoi partitions). 
In lfT6l - lfT8l . distributed algorithms to calculate the standard Voronoi partitions are provided. 
For example, in [[T6ll . explicit stopping criteria are proposed for a generator to calculate its 
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own Voronoi cell without message broadcasting or routing. In contrast, the work in [[TTII . IfTSl 
requires explicitly broadcasting generators or geographic routing. In our work, we assume that 
generator information is available (via either direct sensing or communication). As discussed in 
Section IIII-Bl the energy-based Voronoi partition can also be obtained via a slightly modified 
metric, which happens to be an additively weighted metric as in 11191 . Methods such as Il20ll , 11211 
can then be used to calculate the Voronoi partition based on the modified metric in a centralized 
fashion; however, to the best of our knowledge, there is no known distributed algorithm on 
calculating such partitions. 

The contributions of this work are the following, i) An energy metric is proposed to study 
Voronoi partitions, which arises naturally in battery powered mobile vehicle applications in flow 
environments. To the best of our knowledge, this is the first work on Voronoi partition based on 
the minimum energy required for a vehicle to move from one point to another. In contrast, the 
traveling time based metric has been studied extensively for constant flows [Hl-Q, piecewise 
constant flows and time varying flows [91]. ii) In addition to deriving the lower bound on 
the set of Voronoi neighbors using a disk-based lower approximation of Voronoi cells, an upper 
bound on the set of Voronoi neighbors is proposed utilizing asymptote-based lower and upper 
approximations of Voronoi cells. When deriving the upper bound, we introduce a dominance 
relation among Voronoi generators and provide a complete characterization for the dominance 
relation, iii) Since the upper bound is essential for a vehicle to compute its own Voronoi cell, 
we propose an efficient algorithm based on sorting generators (i.e.. Algorithm ^ besides the 
method based on checking generators sequentially (i.e.. Algorithm [T]). iv) To handle dynamically 
moving vehicles (namely, dynamically changing generators), we introduce a dominance graph 
for recomputing Voronoi cells only when absolutely necessary, which potentially avoids the 
recomputation due to any single change of the set of generators. 

The paper is organized as follows. In Section |IIl we define the minimum energy metric and 
formulate the Voronoi partition problem. Then we characterize the minimum energy metric and 
study the Voronoi partition for two generators in Section [nil In Section |IVl we propose a disk 
based approximation for Voronoi cells and provide a lower bound on the set of Voronoi neighbors. 
To facilitate the calculation of Voronoi partitions in a distributed fashion, we study an asymptote 
based approximation of Voronoi cells, introduce the dominance relation and its characterization, 
and provide an upper bound on the set of Voronoi neighbors in Section |Vl In Section |VI] we 
propose algorithms to calculate the upper bound, and show that the average number of generators 
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in the upper bound is very small via simulations in Section IVIll Finally, we summarize the work 
in Section IVIIIl 

II. Problem Formulation 

In the Cartesian coordinate system, the studied flow environment is described by R^. The 
constant velocity field is a mapping v : (x y)'^ E H"^ ^ {B 0)^, where i? is a positive constant. 
A vehicle runs at speed U = {U^ Uy)^ relative to the velocity field, and then the dynamic of 
the vehicle in the flow environment can be described by 

^ = U.^ + B, (1) 

^ = U (2) 
dt ■ ^ ^ 

We assume that vehicles can run against the flow. 

To study Voronoi partitions, we introduce the following (pseudo)-metric. 

Definition 1 Given two points and in the flow environment R^, the energy metric J{p^,p'^) 
is defined as J{p^,p^) = min jj^ U^Udt, where tf is free, U satisfies Eqs. ([U) and (O, and 
x(0) = Xpi (i.e., the x coordinate of p^), y{0) = ypi (i.e., the y coordinate of p^), x{tf) = 

Xp^, y{tf) = 2/p2- 

The energy metric J{p^,p'^) is the minimum amount of energy required for the vehicle to 
move from its initial location p^ to its final location among all possible controls. Note that 
there is no explicit constraint on U ; however, as shown in Remark [H the optimal control U that 
achieves J{p^,p^) is bounded by two times the flow velocity. The explicit expression for the 
energy metric is derived in Section Hill With this energy metric, now we can define the following 
Voronoi partition. 

Definition 2 Let P = {p^.p'^, ■■■,P^} C R^ be a set of distinct points, where n > 2. We call 
the region given by 

V{p') = {peR'\ J{p\p) < J{p>,p) for J ^ z, J G In] 

the energy-based Voronoi cell associated with p*, where /„ := {1, 2, n}, and the set given by 
V = {V{p^), ^(p^), V{p^)} the energy-based Voronoi partition generated by P. 
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Note that if J{p\p) is replaced with dpip := a/ (xpi — XpY + {upi — VpY, the partition is called 
a standard Voronoi partition. For simplicity, we use Voronoi partition to refer to the energy-based 
Voronoi partition in the rest of the paper. 

Besides calculating the Voronoi partition V given the set of points P, we are especially 
interested in calculating each Voronoi cell V{p'^) for i = l,...,n given F = P \ {]?*}. For 
example, if is the location of a vehicle in the constant flow environment, V{p'^) can be 
interpreted as the set of points that can be reached by with fewer energy consumption than 
by any other vehicle. If a task (e.g., taking measurements) has to be done at a point p belonging 
to V{p^) and vehicle is assigned to the task, the energy consumption is minimized for this 
task. In this context, the challenge of computing the Voronoi cells lies in the fact that since the 
vehicles are moving, the generators p^ for i = l,2,...,n are constantly changing. Without loss 
of generality, we formulate the following Voronoi cell calculation problem. 

Problem 1 Given a fixed point p^ and a set of points P = {p'^,p^, calculate the Voronoi 

cell V{p^) as defined in Definition |2l 

III. Energy-Based Metric and Voronoi Partition with Two Generators 

In this section, we first study the minimum energy control problem and provide an expression 
for the metric J{p^,p^), and then derive the Voronoi boundary between two generators. 

A. Energy Metric: Expression 

The energy metric in Definition [T] is given below. 

Proposition 1 Given two points p^ and p^ in the flow environment with the velocity field v 

satisfying Vx{x,y) = B > and Vy{x,y) = 0, the minimum energy J{p^,p^) = min jj^ U'^Udt 
is 

J{p^,p^) = 2B{dpip2 + Xpi - Xp2) , (3) 
and the optimal control is U{t) = —\ 

C, = 25(l + ^^^^), = ^^^i^^i^, t, = ^. (4) 
dpip2 dpip2 B 

Proof: Refer to the Appendix. ■ 



Ci 

for t G [0, tfl, where 

C2 
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Remark 1 Note that the quantity Jij^.p^) is not a real metric because i) J{p^,p^) = does 
not imply = p^, and ii) J{p^,p'^) is not the same as J(p'^,p^) in general. More specifically, 
if Upi = yp2 and Xpi < Xp2, Eq. ^ reduces to 2B(xp2 — Xpi + Xpi — Xp2) = 0. This is consistent 
with the fact that no control is necessary if lies downstream of p^. In addition, it can be 
verified that only when Xpi = Xp2, J{p^,p^) = J{p^,p^). The magnitude of U satisfies \\U\\ = 
= B. 2+ '^^y'p^^ < 2B. Therefore, the optimal control is bounded. 

B. Voronoi Partition: Two Generators 

Based on the metric J{p^ ,p^) given in Eq. ([3]), we now study the Voronoi partition given two 
generators p^.p^. Following Definition [21 we have V{p^) = {p E | J{p^,p) < J{p'^,p)}- 
Using Eq. ([3]), J{p^,p) < J{p^,p) can be rewritten as 

2B{dpip + Xpi — Xp) < 2B{dp2p + Xp2 — Xp) , 

dp^p ~\~ Xpi ^ dp^p H~ Xp2 ^ (5} 
dp^p dp2p ^ Xp2 Xpi , (6) 

where we obtain Eq. ([5]) because B > 

Based on Eq. Q, we can also use the metric dpip + Xpi to obtain the same Voronoi partition. 
In [fT9ll , this metric falls into the category of additively weighted distances. Therefore, given 
the set of points P, the Voronoi partition in Definition [2] can be calculated using existing 
methods such as EOl . ETI . Straightforward application of such methods to Problem [T] can 
be very inefficient because all Voronoi cells have to be computed in order to just obtain V{p^). 
Another simple idea to solve Problem [T] is that we consider one point at a time and keep refining 
V{p^) until all points have been taken into account. However, as we show in Section |Vl only 
a subset of points are necessary for calculating V{p^). More details on comparing different 
methods to solve Problem [T] are provided in Section IVIIi 

Now we can rewrite V{p^) as V{p^) = {p E^ \ dpip — dp2p < Xp2 — Xpi}, and the boundary 
between p^ and is B{p^,p^) := {j» G | dpip — dp2p = Xp2 —Xpi}. Without loss of generality, 
we assume that Xpi < Xp2 and ypi < yp2. Depending on the relative position between p^ and p^, 
there are three different cases for the Voronoi cells and the boundary between these cells. 

Case I: Xpi < Xp2 and ypi < yp2. For any point p on the boundary we have dpip — dp2p = 
Xp2 — Xpi > 0. It can be verified that the boundary is a hyperbolic curve. To derive an equation. 
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P 



a 



X 




(a) Original coordinates. 



(b) Transformed coordinates. 



Fig. 1. Voronoi partition generated by {p^,p^} satisfying Xpi < Xj,2 and i/pi < yp2. The red dashed hne is the boundary. 



we first transform the coordinate from (x, y) to (x', y') sucli tliat tlie origin is at ( ^ ^ ^ , ' ^ ) 
and the positive x' direction is from to p^. Essentially the transformation involves shifting the 
origin and rotating the x,y axes. As shown in Fig. [Ha), we use a to denote the rotation angle 
Zp^p^p"^, where p^p'^ is parallel to the x axis, and we have tana = ^pL^pL_ Then for any point 
p with coordinates {xp,yp), its coordinate in the {x',y') plane is given as 



x„ = [Xp ) cos a + [Vp , 



sma 



y'p 



1 Xri 



Xpi + Xp2 , y„i -\- 

) sm a + [yp — z — — ) cos a 



(7) 
(8) 



2 . ■ 2 
In the transformed coordinate, the boundary is shown as the red dotted hyperbola in Fig. \V[h), 

and can be described by the equation -^^^ = 1, where a = ^^^^-j^, c = b = 

Vc^ - a2 = V_V ^ and x'p > a. Note that in Fig. \V[h), the coordinate for p* (namely, the 
intersection point between the boundary and the x' axis) is (a, 0). Therefore, in the original 
coordinate, the boundary can be described as 



X 1 -\-x 9 y 1 -\-y 9 X 1 -\-x n y -i -\-y o 

{{xp- P\ P' )cosQ+(yp- P\ P' )sina)2 (-(xp- P\ p' )sina+{yp- ^\ ^' )cosaf 



X n —X 1 
P P \2 



with the constraint that 



•jC I OC ip^: 



) cos a + {yp 



sm a > 



2 / v-/- 2 ' - 2 

Now we apply the above equations of the boundary to specific scenarios. 
Case II: Xpi = Xp2 and ypi < yp2. In this case, a = |, Eq. ^ reduces to yp 



1 (9) 



(10) 



and Eq. (flOl) holds trivially. In other words, the boundary is a perpendicular bisector of the line 
segment p^p"^, as shown in Fig. [2 a). 
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(a) Xpi 


= Xp2 and Hpi < yp2 




(b) 


Xpi < Xp2 and j/^i = yp2. 





Fig. 2. Voronoi partition generated by {p^,p^}. The red dashed line is the boundary. 

Case III: Xpi < Xp2 and ypi = yp2. In this case, a = 0, Eq. ^ reduces to yp = ^^''^ = yp2, 
and Eq. (flOl) reduces to Xp > Xp2. The boundary is a half line given by {p G | > Xp2, yp = 
yp2}, as shown in Fig. I^b). 

It is straightforward to obtain similar equations for the cases with Xpi > Xp2 and/or ypi > yp2. 

IV. Disk-Based Lower Approximation of Voronoi Cells and Lower Bound on 

Voronoi Neighbors 

In this section, we first study disk-based lower approximation of Voronoi cells, and then derive 
a lower bound on the set of Voronoi neighbors. 

A. Disk-Based Lower Approximation 

The disk-based lower approximation of V{p'^) is given as D(p\rpi) = {p G R.^ | dpip < r^^}, 
i.e., a disk centered at with radius r-pi, such that D(p\rpi) C V{p^). 

We first study the case with two points p^ and satisfying Xpi < Xp2. Since in general the 
boundary between the two Voronoi cells is a hyperbola as shown in Fig. [Hb), the radius for 
p^ can be chosen to be rpi = c + a = + ^ j^j^g radius for p'^ can be chosen 

to be rp2 = c — a = — I — ^ ^ ^ . Note that the hyperbola boundary intersects with the 
disk D{p^,rp2) only at one point (namely, the point p* in Fig. [IJb); we denote the point as 
because the focus of the hyperbola is the same as the center of the disk and the 
eccentricity of a hyperbola is larger than \. We have rp2 = dp2p*i^pi p2y It can be verified that 
the hyperbola boundary intersects with the disk D(p^,rpi) also only at the point p*(p^,p'^) and 
Tpi = dpip*(^pi^p2y In the original coordinate, we have p*{p^,p'^) = p^ + ^{P^ — P^)^ where 

X 2—2; 1 , d, I 2 

a = " and c = 
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If Xpi = Xp2 and ypi < yp2 (i.e., Case II as discussed in Section IIII-BI) . then we have 
'^pi = ^p'^ = If Xpi < Xp2 and i/pi = yp2 (i.e., Case III as discussed in Section ITlI-BI) . then 
we have Vpi = dpip2 and rp2 = 0. It can be verified that the above results also hold if Xpi > Xp2. 

If there are n points we can choose the radius Vpi = minjg{i^2,...,n}\{i}(-^ + 

^^^^2^). Therefore, D{p\rp^) CV{p'). 

B. Lower Bound on Voronoi Neighbors 
Now we focus on V{p^) and let 

^nip') = argmin (%i + ^Zli^) . (H) 

In general, Moip^) could be a set with multiple elements. Before analyzing Md{p^), we first 
define the set of Voronoi neighbors of p^. 

Definition 3 Given a fixed point p^ and a set of points P = {p'^,p^, a point p G P is a 

Voronoi neighbor of p^ if V{p^) fl V{p) is non-empty and non-trivial (i.e., not a single point). 
We use A/V(p^) to denote the set of Voronoi neighbors of p^. 

Remark 2 Note that in our setting, V (p^) CiV (p) could be empty, a single point, a line segment, 
or part of a hyperbola. In our definition of Voronoi neighbors, we treat two points as neighbors 
only when the intersection of their Voronoi cells is nonempty and nontrivial (i.e., not a single 
point). This is because we are primarily interested in calculating Voronoi cells and ruling out 
this trivial case does not affect the calculation. A definition of Voronoi neighbors in the same 
spirit is used in [[9l. <) 

Now we are ready to state the relationship between Moip^) and A/V(p^). 

Theorem 1 Given a fixed point p^ and a set of points P = Md{p^) ^ -^vip^)^ 

where J^d{p^) is defined in Eq. (fTT]) . 

Proof: Refer to the Appendix. ■ 
Now it is clear that Moip^) is a lower bound on the set of Voronoi neighbors of p^. Even 
if a point p'' may not minimize the radius of the disk that lower approximates V{p^), p^ can 
still be a Voronoi neighbor of p^ . Thus Moip^) could potentially be augmented to a larger set 
Md{p^) while still satisfying Moip^) ^ J\fv{p^). It can be verified that the following sufficient 
condition for testing if p'^ is a Voronoi neighbor of p^ holds. 
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Proposition 2 Given a fixed point and a set of points P = a point 

for k e {2,3,...,n} is a Voronoi neighbor of p^ if for any / G {2,3,...,n} \ {k} we have 

J{p'',p*{p^,p^)) < J{p\p*{p^,p^)). 

Proof: Refer to the Appendix. ■ 
Essentially, the condition verifies that for a point p'^, the energy required to reach the specific 
point p*(p^,p'') (which belongs to the boundary between p^ and p^) from any point other than p'' 
is strictly larger than the energy from p'^. Based on this condition, we can construct Md{p^) by 
starting with Nd{p^), and adding a point p^ to Moip^) if the condition in Proposition [2] holds. 

V. Asymptote-Based Approximations of Voronoi Cells and Upper Bound on 

Voronoi Neighbors 

In this section, we first propose asymptote-based lower and upper approximations of Voronoi 
cells, and then introduce a dominance relation to upper bound the set of Voronoi neighbors. 

A. Asymptote-Based Approximations of Voronoi Cells 

Since in general the boundary is a hyperbola, another way to approximate the Voronoi cells 
is to use the asymptotes of the hyperbola. Here, we focus on two points p^ and satisfying 
Xpi < Xp2 and i/pi < yp2. In the x' — y' plane of Fig. [flb), the equation for the asymptote li (or I2) 
is y' = ^x' (or y' = —^x'), where a = ^^p^ and b = hH-liL^ \[ can be shown that the region 
described by D[^^^^{p^\p^) = {{x' y'f G | a;' < ^^y> if y' > 0,x' < -fy' if y' < 0} satisfies 
-^iower(^'^b^) — ^'(P^)' where V'{p^) is the Voronoi cell V{p^) in the transformed x' — y' plane. 
Going back to the original x — y plane, we have -Diower(p^lp^) C V{p^). At the same time, we 
get an upper approximation for V{p'^) as -Dupper(j'^b^) = \ -Diower(j'^b^)- 

To obtain an upper approximation for V{p^), we use (which is parallel to li) and (which is 
parallel to I2) that pass through the point p* in Fig.[Ub). The equation for I3 (or U) is y' = ^{x'—a) 
(or y' = —^{x' — a)). It can be shown that the region described by D'^ppsr{p^\p^) = {{x' y'Y £ 
I ^/ < ay, + a if > 0,a;' < -fy' + a if y' < 0} satisfies V'{p^) C D'„^^^Xp^\p'^) . Going 
back to the original x — y plane, we have V{p^) C D^ipp^^{p^\p'^). At the same time, we get a 
lower approximation for V{p'^) as -Diower(p^) = \ -Dupper(p^)- 

If there are n points, we can lower approximate V{p'') using Mower(p*) = ^j<^in\{i}D\ower{p^\p^'), 
and upper approximate V{p') using Kpper(p*) = ^jeir\{i\D^w<^r{.P'\p'), where D{p^\f) is the 
approximation of V{p^) given the point p^ . 
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Since li, . . . ,1^ play a very important role in the approximations, we study their equations in 
the X — y plane. Here we focus on li and I2 since ^3 (or l^) is parallel to li (or I2). For li, we 
are interested in y' = ^x' for y' > 0, while for I2, we are interested in y' = —^x' for y' < 0. 

Proposition 3 Given two points and p"^ in the flow environment satisfying Xpi < Xp2 and 
Z/pi < yp2, two asymptotes of the hyperbolic boundary between the Voronoi cells of p^ and 
are 

ypi+yp2 Xpi+Xp2 .Vpi+Vp^ 

y — = [X —) tan 2a with y > — — , (12) 

2 2 2 

where a = arctan ^p—^pL and 

y= y^'+y^' anda:> ^^^+^^^ . (13) 
y 2 - 2 

Proof: Refer to the Appendix. ■ 

Remark 3 If Xpi < Xp2 and t/pi = yp2, which implies that a = 0, then Eq. (fT2l) becomes 
?/ = VitV ^ 2; > ^p'^~]^^p^ because x' > 0, while Eq. (fT3l) remains the same. In this case, 
the two asymptotes coincide and form the exact boundary as shown in Fig. I2b). If Xpi = Xp2 
and ypi < yp2, which implies that a = |, then Eq. (fT2l) becomes y = hl^!lpL^ and x < ^''^ ^^''^ 
because y' > 0, while Eq. (fT3l) remains the same. In this case, the two asymptotes form a straight 
line and are also the exact boundary as shown in Fig. [2ta). ^ 

Remark 4 Note that the two asymptotes pass through the middle point of p^p'^. One is always 
parallel to the x axis, while the other has the slope tan2Q;. <C> 

Similarly, we can obtain asymptote equations if Xpi > Xp2 and/or ypi > yp2. In the next 
subsection, we introduce a dominance relation and provide conditions to check the dominance 
relation, in which the asymptote equations prove to be useful. 

B. Dominance Relation and its Characterization 

When we calculate the Voronoi cell V{p^) in Problem [T] by considering point p"^ first and then 
p^, it is possible that the Voronoi cell of p^ is not strictly refined when considering p^ given 
p^; in this scenario, p^ is not necessary to compute V{p^), and potentially the calculation of 
V{p^) can be done more efficiently. Now we introduce a dominance relation which captures this 
scenario. 
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Definition 4 Given a fixed point p^, and two points p^,p^ (that are different from p^), we say 
p^ dominates p^ (denoted as p^ >- p^) if V{p^ \ p^) = V{p^ \ p'^.p^), where V{p^ \ p^) = 
{p e I J{p^,p) < J{p^,p)} and V{p^ \ p'^,p^) = G | J{p^,p) < J{p^ ,p)-, J{p^ ,p) < 
J(p^p)}. 

By definition, p^ dominates itself; if p^ dominates p^, then only p^ matters when p^ calculates 
its Voronoi cell (this is proved more generally in Theorem [3]). Given a fixed point p^, to check 
if p^ dominates p^, there are four scenarios: 

• Scenario A: y^i < yp2. The iff condition is given in Theorem [2l 

• Scenario B: ypi > yp2. The iff condition is given in Corollary [H 

• Scenario C: ypi = yp2 and Xpi > Xp2. The iff condition is given in part (a) of Proposition H 

• Scenario D: ypi = yp2 and Xpi < Xp2. The iff condition is given in part (b) of Proposition IH 

Theorem 2 Given a fixed point p^, and two points p'^,p^ that are different from p^ and satisfy 

ypi < yp2, p2 y p3 iff > ^^2^ and 

{yp3 - ypi) X {xp2 -Xpi) < {Xp3 - Xpi) x {yp2 - ypi) . (14) 

Proof: Refer to the Appendix. ■ 
Similarly, we can prove the following result for the case ypi > yp2. 

Corollary 1 Given a fixed point p^, and two points p'^,p^ that are different from p^ and satisfy 

ypi > yp2, p2 y p3 jff < ^^2, and (?/p3 - Vpi) x {xp2 - Xpi) > (xps - a;pi) x {yp2 - ypi). 

If ypi = yp2 for points p^ and p^, the following result can be verified. 

Proposition 4 Given a fixed point p\ and two points that are different from p^ and satisfy 

ypi = yp2 and 

(a) Xpi > Xp2, p^ >- p^ iff yps ^ yp2, or ypz = yp2 and XpS < Xpi. 

(b) Xpi < Xp2, p^ >- p^ iff ypS = yp2 and Xp3 > Xp2. 

Remark 5 Fig. [3] illustrates the region of point p^ that is dominated by in Scenario A. There 
are three different cases depending on the x coordinates of p^ and p"^. If Xp2 > Xpi (or Xp2 = Xpi, 
Xp2 < Xpi, respectively), any point in the red dotted (or green solid, blue dashed, respectively) 
region in Fig. [3] is dominated by p^. <C> 
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Fig. 3. Three different cases for Scenario A (y^i < 11^2): Xp2 > Xpi, = x^i, and Xp2 < x^i. 

C. Upper Bound on Voronoi Neighbors 

In this subsection, we propose an upper bound on the set of Voronoi neighbors based on the 
dominance relation introduced earlier. We first show that the dominance relation is antisymmetric 
under certain conditions. 

Proposition 5 (Antisymmetry of Dominance) Given a fixed point p^, and two points p^,p^ that 
are different from p^, if y^i ^ yp2, or ypi = yp2 and Xpi < Xp2, then >- p^ and p^ y p^ implies 
that p^ = p^. 

Proof: Refer to the Appendix. ■ 
Given a fixed point p^, and two points p'^,p^ that are different from p^, if ypi = yp2 and 
Xpi > Xp2 (namely, Scenario C in Section IV-BI) . p"^ >- p^ and p^ >- p^ may not imply p^ = p^. 
In fact, as long as ypi = yp2 = ypi, Xpi > Xp2 and Xpi > XpS, we have p"^ y p^ and p^ y p^. 
The reason is that for any p satisfying yp = ypi and Xp < Xpi, V{p^\p) = {{x y)^ G | y = 
ypi,x > Xpi } which does not rely on the exact location of p. In this case, the Voronoi cell of p^ 
is degenerated. To simplify the discussion, we make the following assumption, which guarantees 
that the Voronoi cell of p^ is nonempty. 

Assumption 1 Given a fixed point p^, and a set of points P = satisfying points 

in P U {p^} are distinct, \/p e P, it holds that yp ^ ypi, or yp = ypi and Xp > Xpi. 

Given a fixed point p^ and two points p^.p^, if p^,p'^,p^ satisfy Assumption [TJ Proposition |5] 
guarantees that there are only three cases in terms of the dominance relation between and p^, 
i.e., y p^, p^ y p^, or p^ and p^ do not dominate each other. 
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(a) Point locations. (b) Voronoi cell. 



Fig. 4. A fixed point = (0 0)"^ and a set P of 11 points in (a), and the (bounded) Voronoi cell of the point in (b). 
Afaip^) consists of points with the red color in (a), which are also Voronoi neighbors (for this example). In (b), the same color 
is used for plotting the point i and the Voronoi boundary between point i and p^. 

Definition 5 Given a fixed point and a set of points P = satisfying Assump- 

tion [H we define J^ciP^) as the set of points satisfying that for any p G J^ciP^) there does not 
exist another point p' G P that is different from p and dominates p. 

Now we are ready to state the relationship between Mg{p^) and Mv{p^)- 

Theorem 3 Given a fixed point p^ and a set of points P = {p'^,p^, satisfying Assump- 

tion [H M/(p') ^Afaip'). 

Proof: Refer to the Appendix. ■ 

Example 1 Let p^ = (0 0)^, and generate a set P of 11 points that satisfy Assumption [H in the 
square [—10, 10] x [—10, 10] as shown in Fig. ID^a). It can be verified that Mcip^) consists of 
points 6, 7, 9, 11 (which are highlighted using the red color). By calculating the Voronoi cell of 
the point p^ (the bounded version is shown in Fig. Stb)), the Voronoi neighbors of p^ are points 
6, 7, 9, 11, which are the same as the set of points MciP^) (for this example). ^ 

VI. Calculation of the Upper Bound: Algorithms 

In this section, we first discuss how to calculate the upper bound (on the set of Voronoi 
neighbors) in Theorem [3] when the set of generators is fixed, and then propose algorithms to 
deal with dynamically changing generators. 



November 2, 201 1 



DRAFT 



15 



Algorithm 1 Simple Upper Bound Calculation 

Input: A fixed point and a set of points P = {p^ , •■•7P"} satisfying Assumption [T] 
Output: Ng{p^) 

1: Initialize Ng{p^) = 0; 

2: for i 2,3, . . .n do 
3: Let sign ~ 1; 
4: for j = 2, 3, . . . n do 

5: If j ^ i and p^ >- p\ set sign = and exit the inner for loop; 

6: end for 

7: If sign = 1, set Maip^) = AAg(p^) U {p'}; 

8: end for 

9: Output A/bb^)- 



A. Algorithms for Calculating the Upper Bound: Static Case 

Since the set of Voronoi neighbors is important for solving Problem [U it is necessary to 
develop algorithms to calculate the upper bound in Theorem [3l One simple algorithm is given 
in Algorithm [TJ The algorithm just checks the condition in Theorem |3l More specifically, the 
variable sign indicates whether the point belongs to Maip^)'- 1 if it does and otherwise. 
Steps 4-6 verify if p' is dominated by some point p': if it is, then the algorithm exits the inner 
for loop and does not add p" into Ng{p^)\ otherwise (i.e., sign is never set to be 0), the algorithm 
adds to Mcip^)- It can be verified that the algorithm has complexity (9(n^). The complexity 
(9(n^) is tight for the scenario in which none of p'^,p^, ...,p'^ dominates any other point, i.e., 
Moip') = P. 

One natural question to ask is whether there exists more efficient algorithms to calculate the 
upper bound. Note that the conditions in Theorem |2] require that t/pS > yp2. Therefore, if we first 
sort the points according to y coordinates, potentially we can obtain a faster algorithm. 

Given the set of points P = {p'^,p^, we first divide the set of points into three groups: 

points of which the y coordinate is larger than ypi (denoted as P+), points of which the y 
coordinate is the same as ypi (denoted as Pq), and points of which the y coordinate is smaller 
than ypi (denoted as P_). 

We first focus on points in P+. If P_|_ is not empty, we sort the points in P+ according to 
the ascending order of y coordinates and according to the ascending order of x coordinates for 



November 2, 2011 



DRAFT 



16 



points that have the same y coordinate. Suppose there are m points in P+ and the sorted point 
sequence is pi,P2, ■■■,Pm- We use the point anchor to track the angle formed by the positive x 
axis and the ray from to the point (if the point is p, we use Zp to denote this angle). Since pi 
has the smallest y, it must belong to Mcip^) due to Theorem [2l and anchor is initialized as pi. 
Now we consider p2: if P2 is not dominated by the point anchor, then add p2 into J^aiP^) ^nd 
update anchor with p2', otherwise, do nothing. We repeat this procedure for points ^3,^4, ■■■,Pm- 

For points in P_, we can simply change the sign of the y coordinates so that we can use the 
procedure for P+ due to Corollary [T] For points in Pq, we add the point with the smallest x 
coordinate into Maip^)- The detailed algorithm is given in Algorithm [21 

The correctness of the algorithm can be proved as below. The reason why P can be divided 
into three subsets P+, Pq, and P_ is that for any point p in each subset, it can only potentially 
be dominated by points in that subset due to the dominance characterizations in Theorem [2l 
Corollary [H and part (b) of Proposition |4l For points in P+, we obtain pi,p2, ■■■,Pm after 
sorting. Due to the condition in Theorem [2l a point pi can only be dominated by points that 
have indices smaller than i. Thus, pi must belong to Ng{p^) because it has the smallest index, 
and we initialize anchor with pi. When considering the point pi with i = 2, 3, m, the variable 
anchor keeps track of the point p that has formed the largest angle Z.p so far. If pi is dominated 
by the point anchor, it cannot belong to Mcip^)- If Pi is not dominated by the point anchor, 
then Pi cannot be dominated by any point in P+ (and we add pi into A/g(p^) and update anchor 
with Pi); this is because 

i) Pi can only be potentially dominated by points that have indices smaller than i, 

ii) for any point p with its index smaller than i, we have yp. > yp so we only need to check 
the second condition in Theorem [2l however, the second condition essentially just compares 
the angle Zp with the angle Zpi. If there exists a point pf^ such that Zp^ > Zpi, then pi 
is dominated by pk- Since the variable anchor keeps track of the point which forms the 
largest angle so far, it is equivalent to compare the angle Z anchor with Zpi, which is the 
same as checking if pi is dominated by the point anchor. 

Because the condition in Corollary \T\ is symmetric to the condition in Theorem [2l we can just 
change the sign of the y coordinates for points in P_ and apply the procedure for points in P+. 
For points in Pq, there is only one point (namely, the one with the smallest x coordinate) that 
cannot be dominated by any other point due to part (b) of Proposition IH 

It can be verified that the algorithm has complexity 0{n\ogn) because the best sorting 



November 2, 2011 



DRAFT 



17 



Algorithm 2 Efficient Upper Bound Calculation 

Input: A fixed point and a set of points P = {p^ , •■•7P"} satisfying Assumption [T] 
Output: Ng{p^) 

1: Initialize Ng{p^) = 0; 

2: Divide the set of points in P into three subsets P+, Pq, and P_; 

3: if |P+| > then 

4: Sort the points in P+ into pi,p2, ■■■,Pm (with m = |P+|) according to the ascending order of y coordinates 

and according to the ascending order of x coordinates for points that have the same y coordinate; 
5: Add pi into Afoip^) and set anchor — pi\ 
6: for i = 2, 3, . . . m do 

7: If Pi is not dominated by the point anchor, add pi into Afaip^) and set anchor = pi\ 

8: end for 

9: end if 

10: If |Po| > 0, add the point with the smallest x coordinate in Pq into Mg{p^)', 

11: If |P_| > 0, let P'_ — {p I [xp — yp)'^ e I"-}' and apply the procedure in Steps 3-9 by replacing P+ with 

P'_ and adding {xp. — yp^)'^ into Ng{p^)', 

12: Output 7Vg(p^). 



algorithms (e.g., merge sort) have complexity 0{n\ogn) and finding points in MciP^) using 
the sorted point sequence has complexity 0{n). 

B. Algorithms for Calculating the Upper Bound: Dynamic Case 

In this subsection, we study how to calculate the upper bound if the set of generators P 
constantly changes. If we apply the algorithms in Section IVI-AI to changing generators, the 
calculation has to be done for any single change in the generator locations. Here, we introduce 
a dominance graph, which builds up on the dominance relation and can be used to calculate the 
upper bound more efficiently. We first look at another property of the dominance relation. 

Proposition 6 (Transitivity of Dominance) Given a fixed point p\ and three points p'^,p^,p^ 
that are different from p^, if p"^ >- p^ and p^ )~ p^, then p^ )~ p^. 

Proof: Refer to the Appendix. ■ 
Under Assumption [H the dominance relation is a partial order since it is 
• Reflexive because \Jp , p )~ p\ 
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• Antisymmetric because >- and >- p' imply that = p^ as shown in Proposition [51 

• Transitive because p" >- p^ and p^ >- p^ imply that p*^ as shown in Proposition [6l 
Based on this dominance partial order, the set of points P induces a directed graph, which we 
call a dominance graph. 

Definition 6 Given a fixed point p^ and a set of points P = satisfying Assump- 

tion [H we define the dominance graph induced by P as G(p^, P) = {V E), where i) V = P is 
the set of vertices, and ii) E is the set of directed edges and there is a directed edge e E E from 
p* E V to p' E V if p^ )~ pP . In addition, p" is called a parent of p^ , and p' is called a child of 
p". 

Since a point can dominate (or be dominated by) multiple points in P, there could be multiple 
output (or input) edges for this point in the dominance graph. However, there is no cycle in 
dominance graphs as shown in the following proposition, which can be proved via contradiction. 

Proposition 7 Given a fixed point p^ and a set of points P = satisfying Assump- 

tion [B the dominance graph G(p^,P) is acyclic. 

Proof: Refer to the Appendix. ■ 
Therefore, dominance graphs are directed acyclic graphs. Given a dominance graph, we are 
interested in finding points that are not dominated by any other point. 

Definition 7 Given a fixed point p^ and a set of points P = satisfying Assump- 

tion [H a point p G P is a neighbor of p^ in the dominance graph (G(p\P) if there does not 
exist another point p' G P that is different from p and dominates p. 

Based on Definition [51 J^oiP^) is exactly the set of neighbors of p^ in the dominance graph 
G(p^, P). The importance of neighbors in the dominance graph is that only neighbors of p^ matter 
when p^ calculates its own Voronoi cell given the set of points P, as shown in Theorem [3l It 
can be verified that the following result holds. 

Proposition 8 Given a fixed point p^ and a set of points P = {p^,p^, ...,p"} satisfying Assump- 
tion [U p is a neighbor of p^ in the dominance graph iff the in-degree of p in the dominance 
graph G(p\P) is 0. 

Example 2 We still consider the setting in Example [TJ p^ = (0 0)^, and the set of 11 points 
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P is shown in Fig. Ufa). The corresponding dominance graph is shown in Fig. |3a). Based on 
Proposition m points 6, 7, 9, 11 are the neighbors of in the dominance graph. Note that there 
is no cycle in the dominance graph, which is consistent with Proposition U\ ^ 

Now we study how to dynamically maintain the dominance graph when inserting or deleting 
points. We assume that there is an upper bound K on the number of points (naturally, the 
number of mobile vehicles serves as this upper bound K). For each point G P, we have 
the vertex Vpi in the dominance graph, and use the following fields to keep track of the vertex: 
i) ID: i with 2 < i < K; ii) x: the x coordinate of p*; iii) y: the y coordinate of p*; iv) 
Parent, a data array of dimension K. Parent(k) = 1 if there is an edge from Vpk to Vpi, 
otherwise; v) Child: a data array of dimension K. Child{k) = 1 if there is an edge from Vpi 
to Vpk, otherwise; vi) No_of_j)arent: the number of parents of vertex Vpi. In addition, there 
is a data array List_of_neighbor of dimension K for to keep track of its neighbors, i.e., 
List_of_neighbor{k) = 1 if p'^ is a neighbor of p^ and otherwise. 

When inserting a point, the point can affect the child and parent fields of the existing vertices, 
can make a neighbor invalid, and itself can become a new neighbor. The details are given 
in Algorithm [3l The input is a dominance graph G(p^,P) with vertices Vp2, ...,Vpn, a list of 
neighbors, and a new point and the output is the dominance graph G{p^, P U {p^~^^}) and 
the updated list of neighbors. Step 1 initializes the vertex Vpn+i. Steps 2-11 update the fields of 
vertices Vp2,Vp3, Vpn+i. At Step 3, p"~^^ is checked against p\ and there are three cases: 

• dominates Then is a parent of p""*"^, and the number of parents of is 
increased by 1. This corresponds to Step 4; 

• p""*"^ dominates p\ Then is a child of p""*"^, the number of parents of is increased by 
1, and cannot be a neighbor of p^. This corresponds to Step 5. 

• Neither nor p^'^^ dominates. In this case, no change is necessary. 

Steps 7-11 examine whether itself is a neighbor of p^. It can be verified that the insertion 
algorithm has complexity 0{n). Algorithm [3] can start with n = 1 (or n > 2), i.e., when the 
dominance graph is empty. 

When deleting a point, the point can affect the child and parent fields of other vertices, and 
can create new neighbors. The details are given in Algorithm IH The input is a dominance graph 
G{p^, P) with vertices Vp2, Vpn, a list of neighbors, and a point p^ to delete, and the output 
is the dominance graph G(j9^, P \ {p^}) and the updated list of neighbors. Steps 1-5 update the 
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Algorithm 3 Dynamic Point Insertion 

Input: Dominance graph G(p^, P) with each vertex Vp for p E P being represented using the fields ID, x, y, Parent, 
Child, and No_ofj>arent, a global data structure List_of_neighbor, and a new point such that Assumption [T| 
holds for the set of points P U {p"^^} 

Output: Dominance graph G(p^, P U {p"^^}) and updated List_of_neighbor 

1: InitiaUze Vp^+i with ID being n + 1, x being Xp^+i, y being j/pn+i. Parent and C/z/W being vectors of all zeros, 

and No_ofj>arent being 0; 
2: for i = 2,3, . . .n do 

3: Check if dominates or dominates using the results in Theorem |2] Corollary [l] and 

Proposition |4l 

4: If dominates increase V^^+i .No_ofjmrent by 1, and set Vpi .Child{n+1) = 1 and Vp^+i .Parent(i) = 
1; 

5: If dominates p', increase Vpi .No_of_parent by 1, and set Vpi .Parent{n + 1) = 1, List_of_neighbor{i) = 

0, and Vp^+i .Child{i) ^ 1; 
6: end for 

7: if V^n+i .No_of_parent^ then 
8: List_of_neighbor{n 1) = 1; 
9: else 

10: List_of_neighbor{n + 1) = 0. 
11: end if 



fields of vertices Vp2,Vp3, ...,Vpj-i,Vpj+i, ...,Vpn. If is a parent of p-', then remove p' from 
p*'s child list; this is done in Step 2. If is a child of p', then remove from p*'s parent list, 
decrease the number of parents of by 1: if the number of parents of is 0, set j»* to be a 
neighbor of p^; this is done in Step 4. If p^ is a neighbor of p^, Step 6 removes from the list of 
neighbors of p^. It can be verified that the deletion algorithm has complexity 0{n). Algorithmic 
can start with n = 2 {or n > 3), i.e., when there is only one vertex in the dominance graph. 

To obtain the set of neighbors of p^ in the dominance graph at any time, we just need to 
check the nonzero entries in List_of_neighbor. 

Example 3 We still use the point p^ and the set of points P in Example [IJ Now we consider 
inserting points. We first insert point 12 (the relative position of which is shown in Fig. Ob)). 
The dominance graph in Fig.[5ta) becomes the graph in Fig. [2c); here, point 12 is not a neighbor 
of p^ in its dominance graph (since it is dominated by points 2, 7 and 8), which implies that 



November 2, 2011 



DRAFT 



21 



Algorithm 4 Dynamic Point Deletion 

Input: Dominance graph G{p^,F) with each vertex Vp for p E P being represented using the fields ID, x, y, 
Parent, Child, and No_of jmrent, a global data structure List_of_neighbor, and a point p' e P to delete, where 
J e {2, 3, . . . , n} 

Output: Dominance graph G(p^, P \ {p^}) and updated List_of_neighbor 
1: for i = 2, 3, . . . n do 

2: If Vp, .Parent{i) = 1, set Vpi.Child{j) 0; 

3: If Vpj .Child{i) = 1, set Vpi .Parent {j) = and decrease Vpi .No_of_parent by 1; 
4: If Vpi .No_ofjparent~ 0, set List_of_neighbor{i) = 1; 
5: end for 

6: Set List_of_neighbor{i) = 0, and delete the vertex Vpj . 



the Voronoi cell of remains the same following Theorem [3l Next we insert point 13 (the 
relative position of which is shown in Fig. [5tb)). The dominance graph in Fig. [S^c) becomes the 
graph in Fig. Od); here, point 13 becomes a neighbor of in its dominance graph. In this case, 
points 6, 7, 13 are neighbors of p^ in its dominance graph, and can also be verified to be Voronoi 
neighbors. Now we consider deleting points. We first delete point 5 (the relative position of 
which is shown in Fig. Ob)). The dominance graph in Fig. [Sj^d) becomes the graph in Fig. [Sj^e). 
Since the set of neighbors in dominance graph does not change after deleting point 5, the Voronoi 
cell of p^ remains the same as the one after inserting point 13 following Theorem [3l Next we 
delete point 7 (the relative position of which is shown in Fig. [5tb)). The dominance graph in 
Fig. |2e) becomes the graph in Fig. [5tf). Since point 7 is a neighbor of p^ in its dominance 
graph, point 8 becomes a neighbor of p^ in its dominance graph after deleting point 7, as shown 
in Fig. \5lT). In this case, points 6, 8, 13 are neighbors of p^ in its dominance graph, and can also 
be verified to be Voronoi neighbors. Therefore, when a point is inserted or deleted, the Voronoi 
cell could stay the same provided that the set of neighbors in the dominance graph remains the 
same. In contrast, we have to consider all points again after insertion or deletion if without the 
dynamically maintained dominance graph. This becomes particularly useful when the points are 
the positions of mobile vehicles. <0 
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(d) Dominance graph after inserting (e) Dominance graph after deleting (f) Dominance graph after deleting 
point 13, which becomes a neighbor. point 5, which is not a neighbor. point 7, which is a neighbor. 

Fig. 5. Dynamic maintenance of dominance graph. 

VII. Simulations 

In this section, we first run simulations to study the number of neighbors in dominance graphs 
(namely, the cardinality of A/g'(p^)), and then propose methods to solve Problem [TJ 

A. Simulations for the Upper Bound 

Since J^ciP^) is an upper bound on the set of Voronoi neighbors as shown in Theorem |3l it 
is natural to ask the question of how many points there could be in Mcip^). While the exact 
number depends on the specific relative positions of the set of points P, we can study the average 
number of points in MciP^) if the points are generated randomly. 

Let us first fix to be (0 0)^, select the number of points n to generate, and then generate 
each point with uniform distribution in the square [—1, 1] x [—1, 1] independently from all other 
points while satisfying Assumption [T] In other words, the positions of these n points are i.i.d 
(independent and identically distributed). Given the set of generated n points {p^, ...,p'^,p'"~^^}. 
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Njmber of points 



Fig. 6. Average number of points in Nciv^) for randomly generated points over 5000 trials as a function of the number of 
points. 

we can run either Algorithm [T] or Algorithm [2] to calculate Mcip^)- Then the quantity that we 
are interested in is the expected number of points in Ng{p^)- 

Note that if n = 1, we know that any randomly generated point must be in A/g(j5^). Therefore, 
the expected number is 1. For n > 2, we examine it via simulations. We randomly generate n 
points, run 5000 trials, and calculate the average number of points in Afcip^)- The plot of the 
average number of points in Mcip^) as a function of the number of points is given in Fig. [6l As 
can be seen from the figure, the average number of points in Mcip^) (which is an approximation 
of the expected number of points in MaiP^)) always stays below 4.5. This is also confirmed via 
simulations for more points (though less trials are run for the sake of time). More specifically, we 
choose 60, 80, 100, 120, 140, 160 points, and fix the number of trials to be 1000. For example, if 
|P| = 60 and we run 1000 trials, there are 4.34 points in Maip^) on average, and the histogram 
of the number of points in Ng{p^) is plotted in the left upper corner of Fig. |71 Similarly, if 
|P| = 80, 100, 120, 140, 160 and we still run 1000 trials, the average number of points in Mcip^) 
is below 4.5, which is independent of the number of points that we generate. The histograms also 
have similar shapes for all scenarios. The intuition is that, when an additional point is added, the 
probability (that the number of points in Mcip^) increases) decreases as the number of points 
increases because it is more likely that the point is dominated by other points or other points 
dominate this point; therefore, the expected number of points in Ng{p^) will not blow up and 
is very likely to stay around a certain value when the number of points is large enough. 
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Fig. 7. Histograms of thie number of points (and the average number of points) in Ng{p^) for randomly generated points over 
1000 trials. 

B. Methods for Solving Problem [7] 

As discussed in the previous subsection, the number of points in MciP^) is comparatively 
much smaller than the number of points in P especially when the total number of points is 
large. Therefore, it is much more efficient to calculate the Voronoi cell of the point based on 
the points in MaiP^)- In this subsection, we briefly discuss methods for solving Problem [U 

Given a fixed point and a set of points P, one straightforward way to calculate the Voronoi 
cell of p^ is to consider one point p G P at a time and calculate the boundary based on the 
points that have considered so far; we denote this approach as the naive approach. The upper 
bound based approach is that we first calculate J\fcip^), and then calculate the boundary of the 
Voronoi cell only based on the points in Ng{p^)- The difference between these two approaches 
lies in the fact that the naive approach uses all points to directly calculate the Voronoi cell, and 
the upper bound based approach only uses points in MciP^) to do the calculation. Since the 
calculation of the upper bound is much simpler compared with calculating the Voronoi cell of 
p^ (which involves determining the intersection points of hyperbolas), the additional effort to 
preselect the set of points that are necessary to compute the Voronoi cell is well worthy. 

Intuitively, the ratio of the naive approach's running time to that of the dominance graph 
based approach will roughly be the ratio of the total number of points to the number of points 
in Mcip^), which is denoted as IZ, i.e., IZ = jfj^^^p^- To calculate the Voronoi cell, we can 
use any method that is capable of calculating Voronoi cells for additively weighted metrics. For 
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(a) Point locations. 



(b) Voronoi cell. 



Fig. 8. A fixed point = (0 0)^ and a set of 32 points in (a), and the (bounded) Voronoi cell of the point in (b). 



example, the most efficient method is the sweepline algorithm proposed in [|20l : however, the 
algorithm computes the Voronoi cell for every point instead of just one single point, and the 
implementation is complicated since queuing mechanism is necessary. In our simulations, we 
implement a simple (and less efficient) algorithm for calculating the (bounded) Voronoi cell of 
the point p^. The basic idea is that we consider one point p at a time, and determine which part 
of the boundary between p and p^ contributes to the (bounded) Voronoi cell of p^. 

We generate 32 points as shown in Fig. [8ta), and calculate the Voronoi cell. For the naive 
approach, P' = {1, 2, 3, 32}; for the upper bound based approach, P' = {9, 31, 32} which is 
the set of points in Nn{p^). The Voronoi cells are plotted in Fig. [Mb). Note that the upper bound 
based approach run^j 1.55 seconds while the naive approach runs 18.75 seconds. The running 
time ratio is ~ 12.1, while the ratio TZ is ^ ~ 10.7. Note that the two ratios are close as 

1.55 ' S 

expected. By examining the Voronoi cells, they are exactly the same. In Fig. |9l the ratio TZ of 
the number of all points to the average number of points in Mcip^) is plotted by combing the 
data in Fig. [6] and Fig. |7] (note that for the number of points larger than 60, the plot is based 
on interpolation using the data in Fig. |7]). This plot shows that on average the ratio TZ increases 
linearly as a function of the number of points. In other words, the upper bound based approach 
becomes more and more efficient (on average) as the number of points increases. Note that even 
in the worst case, i.e., the largest number of points in Mcip^) given a fixed number of points, 
the ratio TZ is also increasing as can be verified from Fig. Ul 

'Note that the running time for the upper bound based approach does not include the time necessary to compute Afaip^), 
which is ignorable when the total number of points is large. 
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Fig. 9. Plot of the ratio of the number of points to the average number of points in A/'g(p^). 

VIII. Conclusions 

In this paper, we study the Voronoi partition by introducing an energy-based metric in constant 
flow environments. We provide an explicit expression for this energy metric, and use it to derive 
the equation for the Voronoi boundary between two generators. To facilitate the distributed 
calculation of Voronoi cells, we propose a disk-based approximation which leads to a lower 
bound on the set of Voronoi neighbors, and asymptote-based approximations which lead to an 
upper bound on the set of Voronoi neighbors. When deriving the upper bound, we introduce the 
dominance relation and provide a complete characterization. Simulations are run to evaluate the 
upper bound and its effect on calculating Voronoi cells. The results have potential applications 
to any other setting based on additively weighted metrics. 

There are several future directions. First, we would like to study the Voronoi partition when 
an upper bound on the traveling time of vehicles is imposed. This is motivated by applications 
such as search/rescue in which the time taken to reach a point is also very important besides 
saving energy. Second, we would like to generalize the flow field to more general scenarios, 
such as piecewise constant flows (HI, time varying flows [|9l, or even quadratic flows [[T4l . 
Third, we would also like to incorporate the energy necessary for communication ll22l (such as 
transmitting the location information) and sensing into our energy metric. Last, we would like to 
study approximation techniques for Voronoi cells based on metrics other than the energy metric. 
For example, we would like to extend such results to the power metric as studied in [|3l. 
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Appendix 

Proof of Proposition [1] The minimum energy control problem can be formulated as below: 

min / ^ U^Udt 



x{Q) = Xpi, y{0) = ypi ,x{tf) = Xp2, y{tf) = yp2 . 

The objective of the optimization problem is to find a control U which minimizes the total energy. 
Then the Hamiltonian is iJ = U^U + P^{U + N), where P = (Pi Ps)^ and N = {B Of. Using 
the minimum principle [23], we obtain the following coupled ordinary differential equations 
(ODEs) besides Eqs. © and ©: 

dPi 



dt 
dP2 



, (15) 
. (16) 



dt 

Since U is chosen to minimize the Hamiltonian, we have U = —\P- Plugging U into Eqs. © 
and ©, we get the following ODEs: 

From Eq. (fTSl) . we know Pi is a constant, and let Pi(t) = Ci for t G [0,t/]. Similarly, from 
Eq. (fT6l) . we know P2 is also a constant and let P2(t) = C2 for t E [0, t/]. Then using Eqs. (flTl) 
and (fT8l) , we have 



C C t 
x{t) = Xpi + {B - -^)t, and y{t) = ypi ^ . 



Since x{tf) and y{tf) are given, we obtain the following equations 

Xp2 = Xpi + {B-^)tf , (19) 
C t 

Vp^ = — 1^ ■ (20) 

Since tj is, free and there is no cost imposed on the final state in the optimization problem, 
H\tf = 0. Therefore, we have 

+ c'lP = . (21) 
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By solving Eqs. (fT9l) . (|20l) and (|2T]) . we obtain the unique solution in Eq. dH). Then the optimal 



control is U (t) 



Ci 
C2 



for t G [0,t/], and the minimum energy is given in Eq. dS]). 



Proof of Theorem [T] We prove it by contradiction. Suppose p'' G Moip^), but ^ A/'y(p^). 
Then the point p*{p^,p'') must lie outside of V{p^). In other words, there must exist an edge of 
the Voronoi cell V{p^) that intersects with the line segment p^p*{p^ ,p^) at a point p^ (which is 
different from p*{p^,p^)) satisfying d^ipi < c?pip.(pi pfe). Suppose the edge is due to another point 
p^. Since p*{p^,p'^) is the unique point on the boundary between p^ and that is closest to p^, 
we have c?pip.(pipm) < dpipi. Thus, dpi ^p* (^pi ^p^^ < dpipi < dpip*^^pi^pky That is to say, following 
the disk-based approximation, dpi,p.(pi,p™) (i.e., the radius due to the point is strictly smaller 
than c^pijp*(pijpfe) (i.e., the radius due to the point p''). It contradicts with the assumption that 
p'' G J^d{p^)- Therefore, p'' must be a Voronoi neighbor. 



Proof of Proposition |2] Since p*{p^,p^) lies on the boundary between p^ and p^, we have 

J{p^ ,p*{p^ ,p^)) = J(p'^,p*(p-'^,p'^)). Because for any / G {2, 3, ?t,}\{A;} we have J(p'^,p*(p^,p'^)) < 
J{p\p*{p^ ,p^)), we obtain J{p^,p*{p^,p'')) < J{p'^ ,p*{p^ ,p^)) form G {2,3, Based on 

Definition I2I we know that p*(p^,p^) lies in V{p^). Because for any / G {2,3,...,n} \ {A;} 
we have J {p^ , p* (p^ , p'^)) < J{p\p*{p^ ,p^)), p*{p^,p^) must only belong to the edge of V{p^) 
between p^ and p'^, and the edge of V{p^) between p^ and p'^ is not a single point due to the 
strict inequality. Therefore, p^ must be a Voronoi neighbor of p^. 

Proof of Proposition |3] For /i, we are interested in y' = ^x' for y' > 0. Plugging in the 
expressions for x' and y' in Eqs. ^ and ([8]), we obtain 

— [x — —) sm a + [y ) cos a = 

{[x — —) cos a + (y — 1^ ) sin a) tana , 

which simplifies to be ?/ — '^^^ ^^^^ = (x — ^^^^^^^) tan 2a. > is equivalent to — (x — 
^£^)sina + (y - ^^£4^)cosa > 0. Using x - ^^^^ = "'^A^ ^ ^e obtain (y - 
^s^^)!^ > 0. Since a G (0, f ), > reduces to y > V±V _ 
For /2, we are interested in y' = —^x' for y' < 0. Plugging in the expressions for x' and y' 
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in Eqs. (|7]) and ©, we obtain 



— [x — - — —) sm a + [y — - — —) cos a = 

((x — —) cos a + {y — ^ ) sin a) (—tana) 



which can be rewritten as {y — = 0- Since a G (0, |), we have y = . y' < 

+= 
Y 



is equivalent to — (x — iVifsi) sin a < 0. Therefore, x > ^^^iiffii. 



Proof of Theorem H 

If is the same as p"^, p^ >- p^ and the conditions hold trivially. Therefore, in the following 
proof, we assume that p"^ and p^ are different. 

(If part) Depending on the x coordinates of p^ and there are three cases. 

Case I: Xp2 > Xpi. Since yps > yp2 > y^i, we have (xps — Xpi) x {yp2 — ypi) > {yps — ypi) x 
{xp2 — Xpi) > 0, which implies that XpS > Xpi. Therefore, Eq. (fT4l) can be rewritten as 

Vp^ - Vp^ ^ Vp^ - Vp^ ^22) 

Now there are two cases depending on the y coordinate of p^ and p^: 

• If ?/p3 = yp2, Eq. (|22|) implies that x^s > Xp2 since j*^,^'^ are different. For any point p 

in V{p^ I j)^), we have J{p^,p) < J{p^-,p), i.e., 2B{dpip + Xpi — Xp) < 2B{dp2p + Xp2 — 
Xr, . Since dipirp dp3p I dip^pS — dp3p I x pS x p'z (due to the triangular inequality), i.e., 
(ip2p + Xp2 < (ip3p + Xp3, we have 2B{dp2p + Xp2 — Xp) < 2B{dp3p + Xp3 — Xp). Therefore, 
2B{dpip + Xpi — Xp) < 2B{dp'ip + Xp3 — Xp), i.e., J{p^,p) < J{p^p). Therefore, we have 
Vijy" I p^) C I p^,p^). As I p^.p^) C I p2), we have V{p^ \ p^) = 

V{p^ I p^,p^)\ 

• If |/p3 > yp2, Eq. (f22)) implies that Xp3 > Xp2. For any point p m V{p^ \ P^)^ we have 
P ^ -^upper(P^ I P^)' where -Dupper(p^ I P^) is obtained via shifting the vertex of Aipper(p^ I P^) 
to the point p^. In other words, the two lines that consist of the boundary of D'^^^^^{p^ | p^) are 
y = yp2 with x > Xp2 (namely, the half line p^gs in Fig. (TO]), and y — yp2 = [x — Xp2) tan 2a 
with y > yp2 and a being the angle of Zgip^p^ (namely, the half line p^q2 in Fig. [TOl) . 
Now we look at -Diower(p^ I P^)^ the lower approximation for given p^. Recall that 
the two lines that consist of the boundary of -Diower(p^ I P^) are y = ^p^^^p^ with x > 
'^p^'^^p^ (namely, the half line q^q^ in Fig. [TOl where q^ is the middle point of p'^,p^), and 
y — ^p^^^p^ = (x — fiiiV) tan 2/3 with y > '^p^^J'p^ and /3 being the angle of Lq-^p^p^ 
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(namely, the half line q^^q^ in Fig. [101). Since ^ ^ ^ > ?/p2, ^ ^ ^ > Xp2 and /3 < a (due 
to Eq. (l22l)'). we have 2/3 < 2a, and -Dupper(P^ I P^) ^ -Diower(p^ I P^)- In other words, for any 
p e I p^), p e Aower(p^ I P^)- Therefore, 2B{dpip + Xpi -Xp) < 2B{dp2p + Xp2 -Xp) < 
2B{dpZp + XpZ -Xp). Thus, V{p'^ \ p^) C V{p'^ \ p^,p^). Since | p'^,p^) C ^(pi | p^), 
we have | p^) = V{p^ \ p'^,p^). 
Case II: Xp2 = Xpi. Since yp2 > i/pi, Eq. (fT4l) implies that > Xpi = Xp2. There are two 
cases depending on the x coordinates of p^.p^: 

• If Xp3 = Xp2, then we have XpS = Xpi and yps > yp2 since p'^,p^ are different. Li this case, 
the boundary between p^ and p'^ is yi2 = ^ shown in Case II of Section IIII-B[ and 
the boundary between p^ and p^ is ?/i3 = ^^^'^ > ?/i2. Therefore, we have V{p^ \ V^) = 
V{p^ I p'^,p^); 

• If Xp3 > Xp2, then we have Xps > Xpi and ?/p3 > yp2. In this case, the boundary between 

p^ and p^ is ?/i2 = ^^^ii^, and the boundary between p^ and p^ is a hyperbola which is 

y 1 3 

above the lower asymptote yis = ^ 2 ^ — ^12- However, the hyperbola will not touch the 
line yi2. Therefore, we have V{p^ \ p^) = V{p^ \ p^,p^). 
Case III: Xp2 < Xpi. There are three cases depending on the x coordinates of p'^.p^: 

• If Xp3 > Xpi, then Eq. (fT4l) holds because yp3 > yp2 > ypi. The boundary between p^ and 

is a hyperbola which is below the upper asymptote yu = ^Ifiii^, while the boundary 
between p^ and p^ is a hyperbola which is above the lower asymptote 1/13 = ^'"^ ^^^'^ ^ ^12- 
Since the hyperbolas will not touch the asymptotes, we have V{p^ \ p^) = V{p^ \ p^,p^); 

• If Xp3 = Xpi, the boundary between p^ and p^ is a hyperbola which is below the upper 
asymptote yi2 = ^^siiM^ while the boundary between p^ and p'^ is the line yi^ = ^^^'^ > 
yi2- Since the hyperbola will not touch the asymptote y^ due to ypi < yp2, we have 
V{p^ I P^) = vip^ I P^,P^)\ 

• If Xp3 < Xpi, it can be shown that V{p^ \ p^) = V{p^ \ p^,p^) by an argument similar to 
the one used in Case I. 

(Only if part) We prove it via contradiction. Suppose the conditions that ?/p3 > yp2 and (?/p3 — 
i/pi) X (Xp2 — Xpi) < (Xp3 — Xpi) X {yp2 ~ypi) do not hold. Then there are three possibilities: I) y^i < 
Vp^, 2) ypi < ypi < yp2, and 3) ypz > yp2 but {ypi -ypi) x (xp2 -Xpi) > (xpi-Xpi) x {yp2-ypi). 

Case 1): ypz < ypi. In this case, ypS < ypi < yp2. The boundary between p^ and p^ is a 
hyperbola which is below the upper asymptote yrs = ^^^'^'^^^ < y^i. Since yp2 > ypi, the boundary 
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between and is a hyperbola which is above the lower asymptote yi2 = ^ ^"^ ^ ?/p^' ^'^^ 
lies strictly above the boundary between and p^. Therefore, we must have Vijy^ \ C 
I V^)- A contradiction to p^ )~ p^. 
Case 2): y^i < < yp2. Then there are two cases depending on the x coordinates of p^.p^: 

• Xp2 > Xpi. In this case, the boundary between p^ and p"^ lies on/above the asymptote 
2/12 = ^Ifiii^. If Xp3 > Xpi (or Xp3 < Xpi), the boundary between p^ and p^ lies on/above 
(or below) the asymptote yn = < and can be arbitrarily close to the asymptote. 
Therefore, we must have V{p^ \ p^,p^) C V{p^ \ p^). A contradiction. 

• Xp2 < Xpi. In this case, the boundary between p^ and p'^ lies below the asymptote yi2 = 
'^p^'^^p^ and can be arbitrarily close to the asymptote. If XpS > Xpi (or Xpa < Xpi), the 
boundary between p^ and p^ lies on/above (or below) the asymptote yi^ = '^p^^^p'^ < 
and can be arbitrarily close to the asymptote. Therefore, we must have Vijl^ \ p^.p^) C 
V{p^ I P^)- A contradiction. 

Case 3): ypz > yp2 but (i/pS — ypi) x (xp2 — Xpi) > {Xp3 — Xpi) x (yp2 — ypi). There are also 
two cases depending on the x coordinates of p^,p'^: 

• Xp2 > Xpi. The boundary between p^ and lies on/above the asymptote yi2 = ^aii^, and 

can be arbitrarily close to the asymptote. If XpS < Xpi, the boundary between p^ and p^ 

y 1 3 

lies below the asymptote 1/13 = ^ ^ ^ — 2/12 ^^'^ '■^^^ be arbitrarily close to the asymptote. 
Therefore, we must have V{p^ \ p^.p^) C V{p^ \ p^). A contradiction. If XpS > Xpi, 
the upper asymptote of the boundary between p^ and p^ has the slope tan 2/3 with (3 = 
arctan '^^^ > 0. The upper asymptote of the boundary between p^ and has the slope 

p3 pi 

tan 2a with a = arctan '^p^ j^p'^ > q. If {ypi — ypi) x [xp2—Xpi) > (xpS—Xpi) x [yp2—ypi), we 
have (3 > a, which implies that 2/3 > 2a. Therefore, part of the boundary (corresponding to 
the upper asymptote) between p^ and p^ must strictly refine V{p^\p'^), i.e., V{p^ \ p^,p^) C 
V{p^ I P^)- A contradiction. 
. Xp2 < Xpi. Since yp3 > yp2 > ypi, {ypi - ypi) x {xp2 - Xpi) > {xpi - Xpi) x (1/^2 - ypi) 

implies that XpS < Xpi. The lower asymptote of the boundary between p^ and p^ has the 

y 3 ~y 1 

slope tan 2/3 with /3 = arctan p _ p < 0. The lower asymptote of the boundary between 

p3 pi 

p^ and has the slope tan 2a with a = arctan < 0. If {ypi — ypi) x {xp2 — Xpi) > 

(a;p3 — Xpi) X {yp2 — ypi), we have (3 > a, i.e., < — /3 < —a, which implies that 
< —2/3 < —2a. Therefore, part of the boundary (corresponding to the lower asymptote) 
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Fig. 10. Points {p^ ,p^} satisfying x^i < Xp2 and j/^i < -1/^2. The red daslied lines are ttie boundaries. 

between and must strictly refine V{p^\p'^), i.e., V{p^ \ p^,p^) C V{p^ \ p^). A 
contradiction. 

Proof of Proposition |5] If p^ is the same as p^, the result holds trivially. In the following proof, 
we consider the case in which p^ and p^ are different. 

If Upi < yp2 (namely. Scenario A in Section IV-BI ). y p^ implies that yp3 > yp2, and 
(?/p3 —Vpi) X (xp2 — Xpi ) < (xp3 — Xpi ) X (?/p2 —ypi). Since ?/p3 > ?/p2 > ypi , p^ y p"^ implies that 
yp2 > yp3, and {yp2 -ypi)x {xp3-Xpi) < {xp2 -Xpi)x (yps-ypi). From yp^ > yp2 and yp2 > ypi, 
we have yp2 = ypS > ypi. Therefore, (?/p3 — ypi) x {xp2 —Xpi) < {xpS —Xpi) x (yp2 — ypi) implies 
Xp2 < Xpi, and (yp2 — ypi) x (0:^3 — Xpi) < {Xp2 — Xpi) x (?/p3 — ypi) implies XpS < Xp2. Thus, 
Xp2 = Xp3. In summary, = Similarly, we can show that p'^ = p^ for the case ypi > yp2 
(namely. Scenario B in Section IV-BI) . 

If ypi = yp2 and Xpi < Xp2 (namely. Scenario D in Section IV-BI) . p'^ >~ p^ implies that 
Vp^ = and XpS > Xp2. Since ypS = ypi and XpS > Xp2 > Xpi, p^ >- p^ implies that yp2 = ypZ 
and Xp2 > Xp3. Therefore, = p^. 

Proof of Theorem |3] Assumption [T] guarantees that V{p^) is nonempty. Let G J\fv{p^) for 
i E {2, 3, n}, i.e., is a Voronoi neighbor of p^. Then the intersection of the Voronoi cells 
of p^ and is a curve with the set of points being V{p^) fl V{p''). For any p which lies on 
the curve but is not an end point, we have J{p^,p) = J{p\p), and J{p^,p) < Jip' ,p) for any 
3 e {2,3,...,n}\{i}. 

Suppose ^ A/g(p^). Then there exists some p'' for k E {2,3, ...,n} \ {i} such that is 
dominated by p'^. Since p'' >- p\ V{p^\p'') = V{p^\p^,p''). Since for the p chosen previously 
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p G V{p^) C V{p^\p''), we have J{p^,p) < J{p^,p). Since V{p'^\p'^) = V{p^\p'^ ,p^), J{p^,p) < 
J{p^,p) implies that J{p^,p) < J{p\p). Therefore, we must have J{p^,p) < J{p^ ,p) < J{p\p) 
(this can be proved via contradiction). Because J{p^,p) = J{p\p), we have J{p^,p) = J{p'',p). 
This contradicts with J(p\p) < J{p^,p) due to J{p^,p) < J (p^ , p) for any j E {2,3, ...,n}\{i}. 
Therefore, Afv{p^) C Afdp^). 

Proof of Proposition |6] If any pair of points in p'^,p^,p^ is the same, or all three points are the 
same, the result holds trivially. In the following proof, we consider the case in which all three 
points are distinct. 

Depending on the relative position of p^ and p^, there are four scenarios as discussed in 
Section IV-Bl 

Scenario A In this case, ypi < yp2. Since p'^ >- p^ and p^ y p^, we have ypi < yp2 < yps < yp4. 
To show p'^ y p'^, we only need to prove that 

(l/p4 - X -^P^) < i^p'^ -Xpi) X {yp2 - ypi) . (23) 

If Xpi < Xp2, p^ >- p^ implies that Xpi < XpS as argued in the proof of Theorem [2l and 
< ^pL^pL- r»3 y p4 jjjipiies that x„i < x„4, and ^^"^ < ^hL^EeL^ Therefore, we have 
V-v < V-V ^hich implies Eq. ([23]). 

If Xpi = Xp2, Xp3 > Xp2 = Xpi. If Xp3 = Xpi, we must have Xpi > Xpi since p^ >~ p^. Then 
the left hand side of Eq. (l23l) is while the right hand side of Eq. (l23l) is nonnegative because 
Xp4 > Xpi and yp2 > ypi. Therefore, Eq. (|23]) holds. If Xp3 > Xpi, then we have Xpi > Xpi. Then 
Eq. (|23]) holds too. 

If Xpi > Xp2, there are three cases depending on the x coordinates of p^,p^. i) If XpS > Xpi, 
then we have Xp4 > Xpi. Then the left hand side of Eq. (l23l) is negative while the right hand 
side of Eq. ( |23l) is positive because Xp4 > Xpi and yp2 > ypi. Therefore, Eq. ( |23] ) holds, ii) If 
Xp3 = Xpi, we must have Xp4 > Xpi. Then Eq. (|23l) holds, iii) If Xp3 < Xpi and Xp4 > Xpi, then 
Eq. holds. If x„3 < x„i and x„4 < x„i, we have VzV < VzV ^nd VzV < VzV _ 
Therefore, we have < V_J^ which implies Eq. (l23l) holds. 

Scenario B In this case, ypi > yp2. It can be proved in a way similar to the one used for 
Scenario A. 

Scenario C In this case, ypi = yp2 and Xpi > Xp2. If yp2 ^ ypS, then y^s 7^ ypi and the relative 
position of p^ and p^ belongs to either Scenario A or B. However, in both cases, p'^ y p^ implies 
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that yp4 7^ i/pi, i.e., i/pi 7^ yp2. Then >- p^. If ?/p2 = yps and XpS < Xpi, there are two cases: 

• Z/p3 7^ yp4. Then yp4 7^ yp2. Therefore, y p'^. 

• yps = ypi and XpA < Xpi. Then we have yp2 = ypA. Therefore, p^ y p^. 

Scenario D In this case, ypi = yp2 and Xpi < Xp2. We have yps = yp2 and XpS > Xp2 due to 
y p^, and yp4, = yps and Xpi > XpS due to p^ y p^. Therefore, p'^ y p'^. 

Proof of Proposition |7] We show the resuk via contradiction. Suppose there is a cycle >- 
pi+i y y pi-+'j y p« in the dominance graph where all points are in P. Then we have 

pi y pi+'j \yy repeatedly applying the transitivity property (i.e.. Proposition (6]). Since we also 
have p"^^ y p^, p^ = p^^'^ due to the antisymmetry property (i.e.. Proposition [S]). Now the cycle 
becomes p" y p*"*"^ y p*"*"^... y y p\ By repeatedly applying the above analysis, we can 

show that = = ... = = p\ Since the points are assumed to be different from 

each other, a contradiction. 
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